Method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs

ABSTRACT

The method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs estimates porosity and water saturation in a porous medium, such as brine-saturated shale, as is common in carbon capture and storage reservoirs, based upon measured electrical conductivity and seismic P-wave velocity. The estimated porosity and water saturation may be used for monitoring carbon dioxide leakage from a carbon dioxide reservoir to the overlying cap rock of the region. Measured electrical conductivity and seismic P-wave velocity data are used by the present method to estimate the porosity and water saturation in the cap rock. If a decrease in water saturation in the cap rock is found, this indicates that carbon dioxide may be leaking up from the carbon dioxide reservoir. An alert signal is then generated to indicate that there may be a carbon dioxide leak.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to carbon capture and storage (CCS), and particularly to a method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs based on measuring electrical conductivity and seismic P-wave velocity of a cap rock above a carbon dioxide storage reservoir.

2. Description of the Related Art

Carbon capture and storage (CCS) (sometimes also referred to as carbon capture and sequestration) is the process of capturing waste carbon dioxide (CO₂) from large point sources (such as fossil fuel power plants), transporting it to a storage site, and depositing it where it will not enter the atmosphere. The waste CO₂ is typically stored in an underground geological formation. The aim of CCS is to prevent the release of large quantities of CO₂ into the atmosphere from fossil fuel use in power generation and other industries. CCS is a potential means of mitigating the contribution of fossil fuel emissions to global warming and ocean acidification. Although CO₂ has been injected into geological formations for several decades for various purposes, including enhanced oil recovery, the long term storage of CO₂ is a relatively new concept.

One of the most important aspects associated with carbon capture and storage is the ability to detect CO₂ leakage from the storage reservoir at an early stage. Seismic and electromagnetic (EM) methods can be used to detect changes in physical properties, such as porosity and fluid saturation of the subsurface layers due to CO₂ leakage. With regard to EM methods, the electrical conductivity (σ) is sensitive to resistive gaseous or supercritical CO₂ replacing a conductive pore fluid in a porous medium. Thus, electromagnetic surveying can be a significant complementary technique to seismic monitoring. While seismic methods can be used to detect gas with limited quantitative information, electromagnetic methods are robust when it comes to predicting gas saturation.

Seismic measurements can provide information about the elastic properties of the storage medium, such as P-wave velocity (V_(P)), while EM measurements can provide information about the EM properties (namely, electrical conductivity) of the subsurface layers. As CO₂ migrates out of the reservoir, the conductivity and the P-wave velocity of the new host rock will change. The integration of rock properties from seismic and EM data can provide better delineation of the leakage pathway than the physical properties from each method individually. Thus, it would obviously be desirable to provide a technique for the joint inversion of porosity and water saturation from electrical conductivity (σ) and P-wave velocity (V_(P)) measurements.

The electrical properties of rocks depend on the pore geometry, the constituent minerals, the type of fluids and their saturations. Hydrocarbons and CO₂ are considered to be electrically resistive. The electrical resistivity of reservoir rocks is highly sensitive to changes in water saturation. This high sensitivity to water saturation in a reservoir can be exploited by EM techniques, where the response is a function of the rock bulk electrical resistivity.

Archie's law is an empirical relation which is used to describe the electrical resistivity of sedimentary rocks (ρ_(t)) (Ω·m) as a function of water saturation, porosity and pore-fluid resistivity (ρ_(w)) as:

ρ_(t) =aφ ^(−m) S _(w) ^(−n)ρ_(w),   (1)

where m is the cementation factor (between 1.8 and 2.0 for sandstones), n is the saturation exponent (typically around 2), a is the tortuosity factor, φ is the porosity and S_(w) is the water saturation. The formation factor, F, is defined as:

$\begin{matrix} {F = {\frac{a}{\varphi^{m}} = \frac{\rho_{0\;}}{\rho_{w}}}} & (2) \end{matrix}$

where ρ₀ is the resistivity of the rock filled with only water (S_(w)=1). The ratio between ρ_(t) and ρ₀ is known as the resistivity index, I, thus:

$\begin{matrix} {I = {\frac{\rho_{t}}{\rho_{0}}.}} & (3) \end{matrix}$

In terms of conductivity (S/m), Archie's law can be written as:

σ=(1−φ)^(p)σ_(s)+σ_(f)φ^(m),   (4)

where σ_(f) is the conductivity of the fluid filling the pore and σ_(s) is the conductivity of the solid minerals. Each phase has its own connectivity and a representative exponent (m and p). Large exponents occur for low connectivity phases, and small exponents occur for high connectivity phases. The exponent p is given by:

$\begin{matrix} {p = {\frac{\log \left( {1 - \varphi^{m}} \right)}{\log \left( {1 - \varphi} \right)}.}} & (5) \end{matrix}$

The Hermance model is modification of Archie's law that has been developed to study the electrical conductivity of the crust and mantle. The relationship relates the bulk rock resistivity to the pore fluid resistivity and porosity as:

σ=(σ_(f)−σ_(s))Ø^(m)+σ_(s)=(1−φ^(m))σ_(s)+σ_(f)φ^(m).   (6)

If σ_(s)→0, then equation (6) reduces to Archie's law.

The porosity of two-constituent composite rock in terms of the lower and upper Hashin-Shtrikman (HS) bounds on the electrical conductivity is given as:

$\begin{matrix} {\varnothing^{- {HS}} = {\left( \frac{\sigma_{s} - \sigma_{HS}^{-}}{\sigma_{s} - \sigma_{f}} \right) - \left( \frac{\sigma_{f} + {2\sigma_{s}}}{\sigma_{HS}^{-} + {2\sigma_{s}}} \right)}} & (7) \\ {{\varnothing^{+ {HS}} = {\left( \frac{\sigma_{s} - \sigma_{HS}^{+}}{\sigma_{s} - \sigma_{f}} \right) - \left( \frac{3\sigma_{f}}{\sigma_{HS}^{+} + {2\sigma_{s}}} \right)}},} & (8) \end{matrix}$

where σ_(HS) ⁻ and σ_(HS) ⁺ are the lower and upper Hashin-Shtrikman bounds for electrical conductivity, respectively.

A self-similar theory has been developed to study the dielectric response of water-saturated rocks based on a realistic model of the pore space. In order to include the local environmental effects around a grain, a self-similar model is generated by envisioning that each rock grain itself is coated with a skin made of other coated spheres. In the self-similar model, the conductivity of the rock satisfies:

$\begin{matrix} {{\varphi = {\sum\limits_{i = 0}^{N}{p_{i}\left( \frac{\sigma - \sigma_{i}}{{2\sigma} - \sigma_{i}} \right)}}},} & (9) \end{matrix}$

where σ_(i) and p_(i) are the conductivity and the volume fraction of the i-th phase, respectively. A water-wet rock that remains percolating for small values of porosity can be obtained from the assumption that water is the starting host material into which infinitesimal amounts of spheres of matrix and fluids are gradually included. This model is in agreement with Archie's law because it preserves the continuity of the water phase. For two constituents (solid and fluid), the solution is given as:

$\begin{matrix} {{\varphi = {{f\left( {\sigma,\sigma_{s},\sigma_{f}} \right)} = {\left( \frac{\sigma_{s} - \sigma}{\sigma_{s} - \sigma_{f}} \right)\left( \frac{\sigma_{f}}{\sigma} \right)^{W}}}},} & (10) \end{matrix}$

where W=⅓ for spherical inclusions. The EM properties of finely plane-layered media can be obtained by using Backus averaging. Considering a plane-layered medium where each layer is homogeneous, isotropic and thin (compared to the wavelength of the electromagnetic wave), then, if the layers interfaces are parallel to the x-y plane, the properties are independent in the x and y directions. The equivalent medium is transversely isotropic and can be described with two components of the conductivity tensor:

$\begin{matrix} {{\begin{pmatrix} \sigma_{11} \\ \sigma_{33}^{- 1} \end{pmatrix} = {\begin{pmatrix} \sigma_{1} & \sigma_{2} \\ \sigma_{1}^{- 1} & \sigma_{2}^{- 1} \end{pmatrix}\begin{pmatrix} p_{1} \\ p_{2} \end{pmatrix}}},} & (11) \end{matrix}$

where it is assumed that there are two thin layers of proportions p₁ and p₂ and conductivities σ₁ and σ₂ . The proportions p₁ and p₂ can be expressed as:

$\begin{matrix} {{p_{1} = {\frac{\sigma_{11} - \sigma_{2}}{\sigma_{1} - \sigma_{2}} = \frac{\sigma_{11}^{- 1} - \sigma_{2}^{- 1}}{\sigma_{1}^{- 1} - \sigma_{2}^{- 1}}}}{and}{p_{2} = {1 - {p_{1}.}}}} & (12) \end{matrix}$

A formulation that is commonly used to estimate the electromagnetic properties of rocks is known as the complex refractive index method (CRIM). For negligible dielectric permittivity, the CRIM model can be expressed as:

σ=[(1−φ). Σ_(i=1) ^(N) ^(solid) f _(i), (σ_(i))^(1/2)+φ·(S _(w)σ_(w) ^(1/2) +S _(o)σ_(o) ^(1/2) +S _(g)σ_(g) ^(1/2))]²,   (13)

where f_(i) is the fractional volume of the i-th solid component; σ_(i), σ_(w), σ_(o), and σ_(g) are the conductivity of the i-th solid phase, water, oil and gas, respectively; and S_(w), S_(o), and S_(g) are the saturations of water, oil, and gas, respectively. The CRIM model has been found to give good results at high frequencies (above ˜0.5 GHz). The CRIM model has also been used for electromagnetic data at low frequencies (˜200 Hz) to obtain electrical conductivity. The electrical conductivity as a function of saturation, porosity and clay content based on the CRIM model is given by:

σ=[(1−φ)·Cσ _(c) ^(1/2)+φ(1−S _(g))σ_(b) ^(1/2)]²,   (14)

where C is the clay content, σ_(c), and σ_(b) are the clay and brine conductivities, respectively, and S_(g) is the gas saturation, FIGS. 2A and 2B show that the CRIM conductivity-porosity model fits the shale and sandy sections of a well in the Gullfaks field of the North Sea, respectively, better than other models.

The CRIM equation is commonly expressed by the following formula:

√{square root over (σ)}=(1−φ)√{square root over (σ_(m))}+φS _(w)√{square root over (σ_(w))}+φ(1−S _(w))√{square root over (σ_(a))},   (15)

where σ_(m), σ_(w), and σ_(a) are the electrical conductivities of the solid matrix, water, and air in the material, respectively. For a rock with a multi-mineral matrix and/or multi-phase fluid, the CRIM model of equation (13) can be used. σ_(w) is typically calculated as:

σ_(w)=0.15×TDS,   (16)

where TDS indicates the total dissolved solids in parts per thousand (ppt) or g/L. FIG. 3 shows a plot of σ in equation (15) as a function of φ and S_(w) in a typical brine saturated shale. Table 1 lists typical values of rock and fluid properties used to generate plots of the electrical conductivity and P-wave velocity, including the plot of FIG. 3. These values are representative of shales.

TABLE 1 Rock and Fluid Properties for Typical Brine-Saturated Shale Rock/fluid property Value Units Quartz volume (ν_(q)) 0.3 Fraction Clay volume (ν_(c)) 1 − ν_(q) = 0.7 Fraction Rock porosity (φ) 0.1 Fraction Rock water saturation (S_(w)) 1.0 Fraction Conductivity of air (σ_(a)) 0 Siemen/meter (S/m) Conductivity of water (σ_(w)) Equation (13) S/m Conductivity of oil (σ_(o)) 0.0001 S/m Conductivity of quartz (σ_(q)) 0.001 S/m Conductivity of clay (σ_(c)) 0.01 S/m Bulk modulus of quartz (K_(q)) 39 Gigapascal (GPa) Bulk modulus of clay (K_(c)) 20 GPa Bulk modulus of air (K_(a)) 1.01 × 10⁵ Pa Bulk modulus of water (K_(w)) 2.25 GPa Shear modulus of quartz (G_(q)) 40 GPa Shear modulus of clay (G_(c)) 15 GPa Density of air (ρ_(a)) 1 kg/m³ Density of water (ρ_(w)) 1000 kg/m³ Density of quartz (ρ_(q)) 2650 kg/m³ Density of clay (ρ_(c)) 2500 kg/m³ Shale critical porosity (φ_(c)) 0.4 Fraction

Thus, a method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs solving the aforementioned problems is desired.

SUMMARY OF THE INVENTION

The method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs estimates porosity and water saturation in a porous medium, such as brine-saturated shale, as is common in carbon capture and storage (CCS) reservoirs, based upon measured electrical conductivity (σ) and seismic P-wave velocity (V_(p)). The estimated porosity and water saturation may be used for monitoring carbon dioxide (CO₂) leakage from a carbon dioxide reservoir to the overlying cap rock of the region. Controlled source electromagnetic (CSEM) surveying is performed repeatedly at close intervals over the carbon dioxide reservoir to measure the electrical conductivity in the cap rock. A seismic survey is performed right after each CSEM survey in order to measure the P-wave velocity in the cap rock at the same locations where the electrical conductivities were measured. The measured electrical conductivity and seismic P-wave velocity data is then used by the method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs to estimate the porosity and water saturation in the cap rock. If a decrease in water saturation in the cap rock is found corresponding to sequential measurements, this indicates an increase in gas saturation, which indicates that carbon dioxide may be leaking up from the carbon dioxide reservoir. An alert signal is then generated to indicate that there may be a carbon dioxide leak.

The method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs includes the following steps: (a) setting an integer i equal to zero; (b) measuring an electrical conductivity σ and a seismic P-wave velocity V_(p) at a cap rock location above a carbon dioxide storage reservoir; (c) calculating a porosity φ associated with the cap rock as:

$V_{p} = \sqrt{\frac{\begin{matrix} {{\varphi \; K_{a}K_{m}K_{w}} + {\left( {\varphi_{c} - \varphi} \right)\left( {K_{m} + {\frac{4}{3}G_{m}}} \right)}} \\ \left\lbrack {{\varphi_{c}K_{m}{K_{w}\left( {1 - S_{w}} \right)}} + {K_{a}\left( {{\left( {1 - \varphi_{c}} \right)K_{w}} + {\varphi_{c}K_{m}S_{w}}} \right)}} \right\rbrack \end{matrix}}{\begin{matrix} {\varphi_{c}\left\lbrack {{\varphi_{c}K_{m}{K_{w}\left( {1 - S_{w}} \right)}} + {K_{a}\left( {{\left( {1 - \varphi_{c}} \right)K_{w}} + {\varphi_{c}K_{m}S_{w}}} \right)}} \right\rbrack} \\ \left\lbrack {{\left( {1 - S_{w}} \right){\varphi\rho}_{a\;}} + {\left( {1 - \varphi} \right)\rho_{m}} + {\varphi \; S_{w}\rho_{w}}} \right\rbrack \end{matrix}}}$

where a water saturation S_(w) of the cap rock is given by:

${S_{w} = \frac{\sqrt{\sigma} - {\left( {1 - \varphi} \right)\sqrt{\sigma_{m}}} - {\varphi \sqrt{\sigma_{a}}}}{\varphi \left( {\sqrt{\sigma_{w}} - \sqrt{\sigma_{a}}} \right)}},$

where σ_(m) is an electrical conductivity of a solid matrix of the cap rock, σ_(w) is an electrical conductivity of water in the cap rock, and σ_(a) is an electrical conductivity of air in the cap rock, K_(a) is a bulk modulus of the air in the cap rock, K_(m) is a bulk modulus of the solid matrix of the cap rock, K_(w) is a bulk modulus of the water in the cap rock, φ_(c) is a critical porosity of the cap rock, G_(m) is a shear modulus of the solid matrix of the cap rock, ρ_(a) is a density of the air in the cap rock, ρ_(m) is a density of the solid matrix of the cap rock, and ρ_(w) is a density of the water in the cap rock; (d) calculating an i-th value of the water saturation S_(wi) of the cap rock as:

${S_{wi} = \frac{\sqrt{\sigma} - {\left( {1 - \varphi} \right)\sqrt{\sigma_{m}}} - {\varphi \sqrt{\sigma_{a}}}}{\varphi \left( {\sqrt{\sigma_{w}} - \sqrt{\sigma_{a}}} \right)}};$

(e) recording the i-th value of the water saturation S_(wi) of the cap rock in non-transitory, computer-readable memory; and (f) if consecutive i-th values of the water saturation of the cap rock decrease, then generating an alert signal to indicate a carbon dioxide leak in the cap rock, and if the consecutive i-th values of the water saturation of the cap rock do not decrease, then setting i=i+1 and returning to step (b).

These and other features of the present invention will become readily apparent upon further review of the following specification.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram illustrating system components for implementing a method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs according to the present invention.

FIG. 2A is a graph comparing conductivity as a function of porosity modeled by the prior art complex refractive index method (CRIM) against other prior art techniques, including a basic arithmetic model, a harmonic model, a geometric model, the Hermance model, a self-similar model, the Hashin-Shtrikman model, Archie's law, and a modified version of Archie's law (indicated as by the “Glover” curve) for data obtained from a shale section of a well in the Gullfaks field of the North Sea.

FIG. 2B is a graph comparing conductivity as a function of porosity modeled by the prior art complex refractive index method (CRIM) against other prior art techniques, including a basic arithmetic model, a harmonic model, a geometric model, the Hermance model, a self-similar model, the Hashin-Shtrikman model, Archie's law, and a modified version of Archie's law (indicated as by the “Glover” curve) for data obtained from a sandy section of the same well as in FIG. 2A.

FIG. 3 is a three-dimensional plot of electrical conductance as a function of porosity and brine saturation modelled by the prior art complex refractive index method (CRIM) for a typical brine-saturated shale.

FIG. 4 is a three-dimensional plot for P-wave velocity as a function of porosity and brine saturation for a typical brine-saturated shale.

FIG. 5 is a three-dimensional plot of inverted porosity as a function of measured P-wave velocity and electrical conductance generated by the method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs according to the present invention.

FIG. 6 is a three-dimensional plot of inverted brine-saturation as a function of measured P-wave velocity and electrical conductance generated by the method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs according to the present invention.

Similar reference characters denote corresponding features consistently throughout the attached drawings.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

The method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs estimates porosity and water saturation in a porous medium, such as brine-saturated shale commonly found in carbon capture and storage (CCS) reservoirs, based upon measured electrical conductivity (σ) and seismic P-wave velocity (V_(p)). The estimated porosity and water saturation may be used for monitoring carbon dioxide (CO₂) leakage from a carbon dioxide reservoir to the overlying cap rock of the region. Controlled source electromagnetic (CSEM) surveying is performed repeatedly at close intervals over the carbon dioxide reservoir to measure the electrical conductivity in the cap rock. A seismic survey is performed right after each CSEM survey in order to measure the P-wave velocity in the cap rock at the same locations where the electrical conductivities were measured. The measured electrical conductivity and seismic P-wave velocity data is then used by the method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs to estimate the porosity and water saturation in the cap rock. If a decrease in water saturation in the cap rock is found, corresponding to sequential measurements, this indicates an increase in gas saturation, which indicates that carbon dioxide may be leaking up from the carbon dioxide reservoir. An alert signal is then generated to indicate that there may be a carbon dioxide leak.

It should be understood that the calculations of the present method can be performed by any suitable computer system, such as that diagrammatically shown in FIG. 1. Data is entered into the system 100 via any suitable type of user interface 116, and can be stored in memory 112, which can be any suitable type of computer readable and programmable memory and is preferably a non-transitory, computer readable storage medium. Calculations are performed by a processor 114, which can be any suitable type of computer processor, and can be displayed to the user on display 118, which can be any suitable type of computer display. The electrical conductivity σ and seismic P-wave velocity V_(p) are fed to the system by any suitable type of electrical conductivity sensor 120, such as a CSEM system or the like, and by any suitable type of seismic P-wave velocity sensor 122.

The processor 114 can be associated with or incorporated into any suitable type of computing device, for example, a personal computer, a programmable logic controller (PLC), or an application specific integrated circuit (ASIC). The display 118, the processor 114, the memory 112 and any associated computer readable recording media are in communication with one another by any suitable type of data bus, as is well known in the art.

Examples of computer-readable recording media include non-transitory storage media, a magnetic recording apparatus, an optical disk, a magneto-optical disk, and/or a semiconductor memory (for example, RAM, ROM, etc.). Examples of magnetic recording apparatus that can be used in addition to memory 112, or in place of memory 112, include a hard disk device (HDD), a flexible disk (FD), and a magnetic tape (MT). Examples of the optical disk include a DVD (Digital Versatile Disc), a DVD-RAM, a CD-ROM (Compact Disc-Read Only Memory), and a CD-R (Recordable)/RW. It should be understood that non-transitory computer-readable storage media include all computer-readable media, with the sole exception being a transitory, propagating signal.

Although several models have been proposed to explain the relation between V_(P), φ and S_(w) in porous media, the Gassmann model is the most commonly used among these models for low-frequency (less than 100 Hz) seismic studies. The Gassmann model is appropriate for CCS because V_(P) measured from surface seismic reflection surveys has frequencies that fall within this range. In this model, the bulk (K_(s)) and shear (G_(s)) moduli of the saturated material are given by:

$\begin{matrix} {\frac{K_{S}}{K_{m} - K_{s}} = {\frac{K_{d}}{K_{m} - K_{d}} + \frac{K_{f}}{\varphi \left( {K_{m} - K_{f}} \right)}}} & (17) \\ {{G_{s} = G_{d}},} & (18) \end{matrix}$

where K_(m), K_(d) and K_(f) are the bulk moduli of the mineral making up the rock, the dry rock, and the pore fluid, respectively, and G_(d) is the shear modulus of the dry rock. The pore-fluid bulk modulus of an air-water mixture is given by Wood's formula as:

$\begin{matrix} {{\frac{1}{K_{f}} = {\frac{S_{w}}{K_{w}} + \frac{\left( {1 - S_{w}} \right)}{K_{a}}}},} & (19) \end{matrix}$

where K_(w) and K_(a) are the bulk moduli of the pore water and air, respectively. The dry-rock bulk and shear moduli can be approximated by the following linear relations in the porosity range φ≦φ_(c):

$\begin{matrix} {K_{d} = {K_{m}\left( {1 - \frac{\varphi}{\varphi_{c}}} \right)}} & (20) \\ {{G_{d} = {G_{m}\left( {1 - \frac{\varphi}{\varphi_{c}}} \right)}},} & (21) \end{matrix}$

where G_(m) is the shear modulus of the mineral making up the rock and φ_(c) is the critical porosity of the rock. The critical porosity of a porous medium is the porosity at which the medium changes its mechanical behavior from a medium in which the mineral grains are load-bearing to a suspension in which the fluid phase is load-bearing. The typical critical porosity for shales is φ_(c)=0.4. For a multi-mineral rock matrix, K_(m) and G_(m) can be estimated using Hill's average of the Voigt and Reuss bounds as:

$\begin{matrix} {K_{m} = {\frac{1}{2} \times \left\lbrack {{\sum\limits_{i}\; \left( {p_{i} \times K_{i}} \right)} + {\sum\limits_{i}\; \left( \frac{p_{i}}{K_{i}} \right)}} \right\rbrack}} & (22) \\ {{G_{m} = {\frac{1}{2} \times \left\lbrack {{\sum\limits_{i}\; \left( {p_{i} \times G_{i}} \right)} + {\sum\limits_{i}\; \left( \frac{p_{i}}{G_{i}} \right)}} \right\rbrack}},} & (23) \end{matrix}$

where p_(i), K_(i), and G_(i) are the fractional volume, the bulk, and shear moduli of the i-th mineral phase of the rock matrix.

The P-wave velocity in the rock at any porosity or saturation is given by:

$\begin{matrix} {{V_{p} = \sqrt{\frac{K_{s} + {\frac{4}{3}G_{s}}}{\rho_{s}}}},} & (24) \end{matrix}$

where ρ_(s) is the density of the saturated rock, given by:

ρ_(s)=(1−φ)ρ_(m)+φ(1−S _(w))ρ_(a) +φS _(w)ρ_(w),   (25)

where ρ_(m), ρ_(w) and ρ_(a) are the densities of the mineral making up the rock, pore water, and air, respectively. The final form of V_(p) (φ, S_(w)) is given as:

$\begin{matrix} {V_{p} = \sqrt{\frac{\begin{matrix} {{\varphi \; K_{a}K_{m}K_{w}} + {\left( {\varphi_{c} - \varphi} \right){\left( {K_{m} + {\frac{4}{3}G_{m}}} \right)\left\lbrack {{\varphi_{c}K_{m}{K_{w}\left( {1 - S_{w}} \right)}} +} \right.}}} \\ \left. {K_{a}\left( {{\left( {1 - \varphi_{c}} \right)K_{w}} + {\varphi_{c}K_{m}S_{w}}} \right)} \right\rbrack \end{matrix}}{\begin{matrix} {\varphi_{c}\left\lbrack {{\varphi_{c}K_{m}{K_{w}\left( {1 - S_{w}} \right)}} + {K_{a}\left( {{\left( {1 - \varphi_{c}} \right)K_{w}} + {\varphi_{c}K_{m}S_{w}}} \right)}} \right\rbrack} \\ \left\lbrack {{\left( {1 - S_{w}} \right){\varphi\rho}_{a}} + {\left( {1 - \varphi} \right)\rho_{m}} + {\varphi \; S_{w}\rho_{w}}} \right\rbrack \end{matrix}}.}} & (26) \end{matrix}$

FIG. 4 shows a plot of V_(P) in equation (26) as a function of φ and S_(w) in a typical brine-saturated shale. It should be noted that this formulation is valid only from rock porosities from zero to the critical porosity of the rock, since the rock will behave like a suspension beyond this porosity and another formulation is required to estimate V_(P) (e.g., Wood's formula).

The method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs includes the following steps: (a) setting an integer i equal to zero; (b) measuring an electrical conductivity σ and a seismic P-wave velocity V_(p) at a cap rock location above a carbon dioxide storage reservoir; (c) calculating a porosity φ associated with the cap rock as:

$V_{p} = \sqrt{\frac{\begin{matrix} {{\varphi \; K_{a}K_{m}K_{w}} + {\left( {\varphi_{c} - \varphi} \right){\left( {K_{m} + {\frac{4}{3}G_{m}}} \right)\left\lbrack {{\varphi_{c}K_{m}{K_{w}\left( {1 - S_{w}} \right)}} +} \right.}}} \\ \left. {K_{a}\left( {{\left( {1 - \varphi_{c}} \right)K_{w}} + {\varphi_{c}K_{m}S_{w}}} \right)} \right\rbrack \end{matrix}}{\begin{matrix} {\varphi_{c}\left\lbrack {{\varphi_{c}K_{m}{K_{w}\left( {1 - S_{w}} \right)}} + {K_{a}\left( \left( {1 - {\varphi_{c}K_{m}S_{w}}} \right) \right\rbrack}} \right.} \\ \left\lbrack {{\left( {1 - S_{w}} \right){\varphi\rho}_{a}} + {\left( {1 - \varphi} \right)\rho_{m}} + {\varphi \; S_{w}\rho_{w}}} \right\rbrack \end{matrix}}}$

where a water saturation S_(w) of the cap rock is given by:

${S_{w} = \frac{\sqrt{\sigma} - {\left( {1 - \varphi} \right)\sqrt{\sigma_{m}}} - {\varphi \sqrt{\sigma_{a}}}}{\varphi \left( {\sqrt{\sigma_{w}} - \sqrt{\sigma_{a}}} \right)}};$

where σ_(m) is an electrical conductivity of a solid matrix of the cap rock, σ_(w) is an electrical conductivity of water in the cap rock, and σ_(a) is an electrical conductivity of air in the cap rock, K_(a) is a bulk modulus of the air in the cap rock, K_(m) is a bulk modulus of the solid matrix of the cap rock, K_(w) is a bulk modulus of the water in the cap rock, φ_(c) is a critical porosity of the cap rock, G_(m) is a shear modulus of the solid matrix of the cap rock, ρ_(a) is a density of the air in the cap rock, ρ_(m) is a density of the solid matrix of the cap rock, and ρ_(w) is a density of the water in the cap rock; (d) calculating an i-th value of the water saturation S_(wi) of the cap rock as:

${S_{wi} = \frac{\sqrt{\sigma} - {\left( {1 - \varphi} \right)\sqrt{\sigma_{m}}} - {\varphi \sqrt{\sigma_{a}}}}{\varphi \left( {\sqrt{\sigma_{w}} - \sqrt{\sigma_{a}}} \right)}};$

(e) recording the i-th value of the water saturation S_(wi) of the cap rock in non-transitory, computer-readable memory; and (f) if consecutive i-th values of the water saturation of the cap rock decrease, then generating an alert signal to indicate a carbon dioxide leak in the cap rock, and if the consecutive i-th values of the water saturation of the cap rock do not decrease, then setting i=i+1 and returning to step (b).

In step (c), the nonlinear equation in φ results in two real solutions. Checking the validity of the solutions using known models showed that only one of the solutions consistently gave the correct porosities of tested model, which is used as the value φ(σ, V_(P), K_(m), K_(w), K_(a), G_(m), φ_(c), ρ_(m), ρ_(w), ρ_(a), σ_(m), σ_(w), σ_(a)). The other solution was found to always give the wrong porosity values, and is thus neglected.

The calculated inverted φ(σ, V_(P), K_(m), K_(w), K_(a), G_(m), φ_(c), ρ_(m), ρ_(w), ρ_(a), σ_(m), σ_(w), σ_(a)) from step (c) and the inverted S_(wi)(σ, V_(P), K_(m), K_(w), K_(a), G_(m), φ_(c), ρ_(m), ρ_(w), ρ_(a), σ_(m), σ_(w), σ_(a)) from step (d) for the case of a brine-saturated shale using the values in Table 1 for the mineral and pore-fluid properties are shown in FIGS. 5 and 6, respectively. The method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs has been applied to 17 σ-V_(p) well-log measurements in a shale section of a vertical well from the Gullfaks field of the North Sea, which was fully saturated with brine. Table 2 provides a list of these σ-V_(P) measurements.

TABLE 2 σ-V_(P) Measurements of Shale Section Measured V_(P) (km/s) Measured σ (S/m) Lithology 2.90 1.52 Shale 2.80 1.51 Shale 2.79 1.53 Shale 2.79 1.52 Shale 2.70 1.51 Shale 2.70 1.48 Shale 2.63 1.55 Shale 2.70 1.54 Shale 2.61 1.60 Shale 2.59 1.50 Shale 2.60 1.42 Shale 2.40 2.12 Shale 2.41 2.20 Shale 2.38 2.30 Shale 2.41 2.35 Shale 2.42 2.21 Shale 2.50 2.40 Shale

Table 3 shows the well-log porosities measured in this depth interval alongside those estimated by the present method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs, with absolute errors shown in the rightmost column. The absolute errors lie between 0.5% and 25% with a mean of 7.35% and a standard deviation of 7.33%. The points with relatively high error may be attributed to local heterogeneities within this shale interval caused by different matrix and/or fluid properties than those used for the inversion.

TABLE 3 Comparison Between Measured and Estimated Porosities Measured φ Estimated φ Absolute Error (%) 0.221 0.229 3.8 0.222 0.228 2.8 0.224 0.231 2.9 0.241 0.229 4.8 0.242 0.228 5.7 0.254 0.225 11.6 0.261 0.233 10.8 0.272 0.232 14.8 0.274 0.239 12.9 0.282 0.227 19.5 0.292 0.217 25.6 0.300 0.294 1.8 0.301 0.302 0.5 0.31 0.312 0.7 0.312 0.317 1.6 0.313 0.303 3.1 0.315 0.322 2.1

Table 4 shows the water saturation values measured in this depth interval alongside those estimated by the present method, with absolute errors in the rightmost column. The absolute errors lie between 0.002% and 0.004% with a mean of 0.003% and a standard deviation of 0.0005%.

TABLE 4 Comparison Between Measured and Estimated Water Saturations Measured S_(w) Estimated S_(w) Absolute Error (%) 1.0 1.00004 0.004 1.0 1.00004 0.004 1.0 1.00004 0.004 1.0 1.00004 0.004 1.0 1.00003 0.003 1.0 1.00003 0.003 1.0 1.00003 0.003 1.0 1.00003 0.003 1.0 1.00003 0.003 1.0 1.00002 0.002 1.0 1.00002 0.002 1.0 1.00002 0.002 1.0 1.00003 0.003 1.0 1.00003 0.003 1.0 1.00003 0.003 1.0 1.00003 0.003 1.0 1.00003 0.003

With the exception of a few samples with relatively high errors, the present method estimated rock porosity within this formation very well. The high-error porosities are clustered in the middle of this shale section, which might be explained by a layer having matrix and/or fluid properties that are slightly different from those used for the inversion. If this is the case, then the background parameters have to be adjusted accordingly. The water saturation results demonstrate much higher accuracy, although readings with a wider range of water saturations are required to confirm this result.

It is to be understood that the present invention is not limited to the embodiments described above, but encompasses any and all embodiments within the scope of the following claims. 

I claim:
 1. A method of monitoring carbon dioxide leakage in carbon capture and storage reservoirs, comprising the steps of: (a) setting an integer i equal to zero; (b) measuring an electrical conductivity σ and a seismic P-wave velocity V_(p) at a cap rock location above a carbon dioxide storage reservoir; (c) calculating a porosity φ associated with the cap rock as $V_{p} = \sqrt{\frac{\begin{matrix} {{\varphi \; K_{a}K_{m}K_{w}} + {\left( {\varphi_{c} - \varphi} \right){\left( {K_{m} + {\frac{4}{3}G_{m}}} \right)\left\lbrack {{\varphi_{c}K_{m}{K_{w}\left( {1 - S_{w}} \right)}} +} \right.}}} \\ \left. {K_{a}\left( {{\left( {1 - \varphi_{c}} \right)K_{w}} + {\varphi_{c}K_{m}S_{w}}} \right)} \right\rbrack \end{matrix}}{\begin{matrix} {\varphi_{c}\left\lbrack {{\varphi_{c}K_{m}{K_{w}\left( {1 - S_{w}} \right)}} + {K_{a}\left( {{\left( {1 - \varphi_{c}} \right)K_{w}} + {\varphi_{c}K_{m}S_{w}}} \right)}} \right\rbrack} \\ \left\lbrack {{\left( {1 - S_{w}} \right){\varphi\rho}_{a}} + {\left( {1 - \varphi} \right)\rho_{m}} + {\varphi \; S_{w}\rho_{w}}} \right\rbrack \end{matrix}}}$ where a water saturation S_(w) of the cap rock is given by: ${S_{w} = \frac{\sqrt{\sigma} - {\left( {1 - \varphi} \right)\sqrt{\sigma_{m}}} - {\varphi \sqrt{\sigma_{a}}}}{\varphi \left( {\sqrt{\sigma_{w}} - \sqrt{\sigma_{a}}} \right)}},$ wherein σ_(m) is an electrical conductivity of a solid matrix of the cap rock, σ_(w) is an electrical conductivity of water in the cap rock, and σ_(a) is an electrical conductivity of air in the cap rock, K_(a) is a bulk modulus of the air in the cap rock, K_(m) is a bulk modulus of the solid matrix of the cap rock, K_(w) is a bulk modulus of the water in the cap rock, φ_(c) is a critical porosity of the cap rock, G_(m) is a shear modulus of the solid matrix of the cap rock, ρ_(a) is a density of the air in the cap rock, ρ_(m) is a density of the solid matrix of the cap rock, and ρ_(w) is a density of the water in the cap rock; (d) calculating an i-th value of the water saturation S_(wi) of the cap rock as: ${S_{wi} = \frac{\sqrt{\sigma} - {\left( {1 - \varphi} \right)\sqrt{\sigma_{m}}} - {\varphi \sqrt{\sigma_{a}}}}{\varphi \left( {\sqrt{\sigma_{w}} - \sqrt{\sigma_{a}}} \right)}};$ (e) recording the i-th value of the water saturation S_(wi) of the cap rock in non-transitory, computer-readable memory; and (f) if consecutive i-th values of the water saturation of the cap rock decrease, then generating an alert signal to indicate a carbon dioxide leak in the cap rock, and if the consecutive i-th values of the water saturation of the cap rock do not decrease, then setting i=i+1 and returning to step (b).
 2. A computer software product that includes a non-transitory storage medium readable by a processor, the non-transitory storage medium having stored thereon a set of instructions for performing monitoring carbon dioxide leakage in carbon capture and storage reservoirs, the instructions comprising: (a) a first set of instructions which, when loaded into main memory and executed by the processor, causes the processor to set an integer i equal to zero; (b) a second set of instructions which, when loaded into main memory and executed by the processor, causes the processor to measure an electrical conductivity σ and a seismic P-wave velocity V_(p) at a cap rock location above a carbon dioxide storage reservoir; (c) a third set of instructions which, when loaded into main memory and executed by the processor, causes the processor to calculate a porosity φ associated with the cap rock as: $V_{p} = \sqrt{\frac{\begin{matrix} {{\varphi \; K_{a}K_{m}K_{w}} + {\left( {\varphi_{c} - \varphi} \right){\left( {K_{m} + {\frac{4}{3}G_{m}}} \right)\left\lbrack {{\varphi_{c}K_{m}{K_{w}\left( {1 - S_{w}} \right)}} +} \right.}}} \\ \left. {K_{a}\left( {{\left( {1 - \varphi_{c}} \right)K_{w}} + {\varphi_{c}K_{m}S_{w}}} \right)} \right\rbrack \end{matrix}}{\begin{matrix} {\varphi_{c}\left\lbrack {{\varphi_{c}K_{m}{K_{w}\left( {1 - S_{w}} \right)}} + {K_{a}\left( {{\left( {1 - \varphi_{c}} \right)K_{w}} + {\varphi_{c}K_{m}S_{w}}} \right)}} \right\rbrack} \\ \left\lbrack {{\left( {1 - S_{w}} \right)\varphi \; \rho_{a}} + {\left( {1 - \varphi} \right)\rho_{m}} + {\varphi \; S_{w}\rho_{w}}} \right\rbrack \end{matrix}}}$ where a water saturation S_(w) of the cap rock is given by: $S_{w} = \frac{\sqrt{\sigma} - {\left( {1 - \varphi} \right)\sqrt{\sigma_{m}}} - {\varphi \sqrt{\sigma_{a}}}}{\varphi \left( {\sqrt{\sigma_{w}} - \sqrt{\sigma_{a}}} \right)}$ wherein σ_(m) is an electrical conductivity of a solid matrix of the cap rock, σ_(w) is an electrical conductivity of water in the cap rock, and σ_(a) is an electrical conductivity of air in the cap rock, K_(a) is a bulk modulus of the air in the cap rock, K_(m) is a bulk modulus of the solid matrix of the cap rock, K_(w) is a bulk modulus of the water in the cap rock, φ_(c) is a critical porosity of the cap rock, G_(m) is a shear modulus of the solid matrix of the cap rock, ρ_(a) is a density of the air in the cap rock, ρ_(m) is a density of the solid matrix of the cap rock, and ρ_(w) is a density of the water in the cap rock; (d) a fourth set of instructions which, when loaded into main memory and executed by the processor, causes the processor to calculate an i-th value of the water saturation S_(wi) of the cap rock as: ${S_{wi} = \frac{\sqrt{\sigma} - {\left( {1 - \varphi} \right)\sqrt{\sigma_{m}}} - {\varphi \sqrt{\sigma_{a}}}}{\varphi \left( {\sqrt{\sigma_{w}} - \sqrt{\sigma_{a}}} \right)}};$ (e) a fifth set of instructions which, when loaded into main memory and executed by the processor, causes the processor to record the i-th value of the water saturation S_(wi) of the cap rock in non-transitory, computer-readable memory; and (f) a sixth set of instructions which, when loaded into main memory and executed by the processor, causes the processor to generate an alert signal to indicate a carbon dioxide leak in the cap rock if consecutive i-th values of the water saturation of the cap rock decrease, and to set i=i+1 and return to the second set of instructions if the consecutive i-th values of the water saturation of the cap rock do not decrease. 